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ABSTRACT 

Using HST/ AGS slitless grism spectra from the PEARS program, we study the stellar populations 
of morphologically selected early-type galaxies in the GOODS North and South fields. The sample 

extracted from a visual classification of the (v2.0) HST/ ACS images and restricted to redshifts 
z>0.4 - comprises 228 galaxies (iF775W < 24 mag, AB) out to z<1.3 over 320 arcmin 2 , with a median 
rcdshift zm=0.75. This work significantly increases our previous sample from the GRAPES survey in 
the HUDF (18 galaxies over ~11 arcmin 2 ; Pasquali et al. 2006b). The grism data allow us to separate 
the sample into 'red' and 'blue' spectra, with the latter comprising 15% of the total. Three different 
grids of models parameterising the star formation history are used to fit the low-resolution spectra. 
Over the redshift range of the sample - corresponding to a cosmic age between 5 and 10 Gyr - we 
find a strong correlation between stellar mass and average age, whereas the spread of ages (defined 
by the RMS of the distribution) is roughly ~ 1 Gyr and independent of stellar mass. The best-fit 
parameters suggest it is formation epoch and not formation timescale, that best correlates with mass 
in early-type galaxies. This result - along with the recently observed lack of evolution of the number 
density of massive galaxies - motivates the need for a channel of (massive) galaxy formation bypassing 
any phase in the blue cloud, as suggested by the simulations of Dekel et al. (2009). 

Subject headings: galaxies: elliptical and lenticular, cD - galaxies: evolution - galaxies: stellar content 



1. INTRODUCTION 

For nearly a century, we have known that the local 
galaxy population follows a bimodal colour distribution, 
with a majority of star-forming galaxies and a minority 
of early-type systems (Hubble 1926; Strateva et al. 2001; 
Blanton et al. 2003). This distribution (which has been 
observed up to z~l, Bell et al. 2004) is commonly in- 
terpreted as the outcome of different star-formation his- 
tories, with blue galaxies undergoing a much more pro- 
longed star-formation activity than the early-types. This 
does not rule out that early-type galaxies may also ex- 
perience occasional episodes of star formation at z~0. 
However, if they do, they produce an amount (in mass) of 
young stars definitively lower than the blue-peak galax- 
ies (Trager et al. 2000; Ferreras et al. 2006; Rogers et 
al. 2007). Early- type galaxies are known to exhibit a 
tight colour-magnitude relation (up to z<l, Stanford et 
al. 1998) which is the direct result of a metallicity se- 
quence where more massive galaxies are also metal-richer 
(Bower et al. 1992; Kodama & Arimoto 1997; Vazdekis 
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et al. 2001; Bernardi et al. 2003). The small scatter in 
their colours and mass-to-light ratios (see e.g. Kelson 
et al. 2000), and their relatively high [a/Fe] (Nelan et 
al. 2005; Thomas et al. 2005) suggest that early-type 
galaxies formed the bulk of their stars at high redshift 
(z>3 van Dokkum et al. 1998; Thomas et al. 2005) 
and on a relatively short timescale (i.e. ~1 Gyr for 
the more massive galaxies, Thomas et al. 1999). While 
these properties can be well reproduced by the "mono- 
lithic collapse hypothesis" (Eggen et al. 1962; Larson 
1974), this scenario does not fit well within the currently 
accepted ACDM galaxy formation paradigm, in which 
galaxies form through hierarchical buildup. In semi- 
analytical models, the hierarchical growth of dark mat- 
ter structures is complemented by the addition of simple 
prescriptions of the baryonic physics which describe star 
formation, chemical enrichment, stellar and AGN feed- 
back (e.g. Scannapieco et al. 2005; Croton et al. 2006; 
Hopkins et al. 2006; Monaco et al. 2007; Hopkins et 
al. 2008; Somerville et al. 2008). The model predictions 
are in reasonable agreement with the observed proper- 
ties of early-type galaxies (see e.g. Kaviraj et al. 2005; 
De Lucia et al. 2006; Khochfar & Silk 2006a,b), but they 
still can not reproduce the observed high [a/Fe] ratios. 
On the other hand, support to the merging scenario is 
given by the detection, at rest-frame UV wavelengths, of 
recent episodes of star formation in early-type galaxies 
(Ferreras & Silk 2000; Kaviraj et al. 2007; Rogers et 
al. 2007; Schawinski et al. 2007; Kaviraj et al. 2008) 
and by the evidence that less massive early-type galax- 
ies increased their mass by 20 - 40% since z=l (Bell et 
al. 2004; Ferreras et al. 2005; Treu et al. 2005a) 

We note that all the knowledge collected on early-types 
is primarily based on the nearby galaxy population. Ex- 
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tending the study to a large number of early- type galaxies 
at intermediate and high redshifts allows us to trace their 
mass-assembly and star formation histories in great de- 
tail and, consequently, to fine-tune the theory of galaxy 
formation and evolution. Ferreras et al. (2009a,b) de- 
rived the evolution of number density and size (i.e. half 
light radius) for a sample of early-type galaxies between 
z=l and z=0, extracted from the GOODS North and 
South fields (see also Ferreras et al. 2005) and compared 
it with recent surveys such as GOODS/MUSIC (Fontana 
et al. 2006), DEEP2 (Conselice et al. 2007; Trujillo 
et al. 2007), COMB017 Borch et al. (2006), 2dFGRS 
(Croton et al. 2005) and SDSS (Shen et al. 2003). They 
found that massive early-types do not exhibit any signif- 
icant change in comoving number density over the past 8 
billion years. On the contrary, there exists a noticeable 
evolution in size, implying an increase of about a factor 
of 2 for early-type galaxies more massive than lO n M0 
and a factor of ^50% for the less massive ones. This 
result is in agreement with Daddi et al. (2005), Con- 
selice et al. (2007), Trujillo et al. (2006,2007), Scarlata 
etal. (2007), van Dokkum et al. (2008) and van der Wei 
et al. (2008). The latter also found that the a-R e rela- 
tion is different for nearby and distant early-type galax- 
ies, suggesting a significant change in the dynamics of 
these galaxies. Given the exposure times needed for ac- 
curate velocity dispersions in these galaxies, it does not 
come as a surprise that the amount of evolution of the 
velocity dispersion is still controversial (see e.g. Cenarro 
& Trujillo 2009; van Dokkum et al. 2009; Cappellari et 
al. 2009). Comparison with the semi-analytical mod- 
els of Khochfar & Silk (2006a, b) - where galaxy growth 
occurs via dissipative major mergers - suggests that the 
constant comoving number density at high M* is due 
to a balance between the "sink" (i.e. loss due to merg- 
ers of massive galaxies producing more massive galaxies) 
and the "source" terms (i.e. gain from mergers at lower 
mass) in the 0<z <1 range. Ferreras et al. (2009a) also 
showed that the slope of the Kormendy relation in the 
0.4<z<1.4 range is consistent with z=0 values (see also 
Ziegler et al. 1999; Waddington et al. 2002). Only the 
average surface brightness is seen to change, according 
to pure passive evolution of old stellar populations. 

The fact that stellar populations in early-type galaxies 
undergo passive evolution is confirmed by their colours 
and photometric spectral energy distributions (SEDs). 
Morphologically-selected early-type galaxies have been 
so far identified in clusters up to z~1.4 and their colours 
have been proven to be consistent with those of z=0 
early-types through pure passive evolution of their stars 
(see e.g. Dressier et al. 1990; Aragon-Salamanca et 
al. 2003; Rakos & Schombert 1995; van der Marel & 
van Dokkum 2007; Cool et al. 2008; Whiley et al. 2008). 
Ellis et al. (1997), Stanford et al. (1998) and Blakeslee 
et al. (2003) showed that early-type galaxies in clus- 
ters at z~0.5 - 1.2 complete their star formation by z>3 
(see also Ferreras et al. 1999; Gladders et al. 1998; 
Kodama et al. 1998; Nelson et al. 2001; van Dokkum 
et al. 2000). The same can be said for the most mas- 
sive early-type galaxies (M* > 10 n M©) in the field at 
the same redshift, while the less massive ones appear to 
sustain a more prolonged star formation (see e.g. Bell et 
al. 2004; van der Wei et al. 2004, 2005; van Dokkum & 
Ellis 2003; Treu et al. 2005a,b). 



Although photometric colours and SEDs are clearly 
the most affordable observables to obtain in terms of 
telescope time, they do not disentangle between age and 
mctallicity effects as accurately as spectra and line in- 
dices. The measurement of the Mg6 index and of the 
Hs and H 7 strength in early-type galaxies up to z=0.83 
has independently provided z>2 - 3 as the redshift at 
which their stars formed (Ziegler & Bender 1997; Kelson 
et al. 2001). Spectroscopy from the ground has made 
also possible to use a set of Fell, Mgll and Mgl lines in 
the rest-frame UV and identify massive early-type galax- 
ies (M* > lO n M ) up to z~2.15 (Cimatti et al. 2004; 
Glazebrook et al. 2004; McCarthy et al. 2004; Saracco 
etal. 2005). In particular, Cimatti et al. (2004) stacked 
together the rest-frame UV spectra of 4 early-type galax- 
ies with 1.6 < z < 1.9 and compared the average with 
synthetic spectra of simple stellar populations to find a 
mean stellar age of about 1 Gyr for solar metallicity (or 
better 0.5 < age < 1.5 Gyr for 2.5 > Z/Z G > 0.4). Sim- 
ilarly, Cimatti et al. (2008) combined the spectra of 13 
early-type galaxies at (z)=1.6 and determined a mean 
stellar age in the range 0.7 - 2.8 Gyr for Z = 1.5 — 0.2Z Q , 
corresponding to a formation redshift zp >2. Along this 
line, Kriek et al. (2008) averaged the rest-frame opti- 
cal spectra of 16 early-type galaxies at 2<z<3 to find 
the typical SED of a poststarburst galaxy with an age of 
about 1 Gyr for solar metallicity. They concluded that 
early-type galaxies are likely to form between zp=2 and 
3 and are not yet in place at z=3 (see also Dunlop et 
al. 1996; Spinrad et al. 1997; Driver et al. 1998; Bram- 
mcr & van Dokkum 2007; Kodama et al. 2007). 

At redshifts z>l, though, the performances of ground- 
based spectrographs is severely limited by the atmo- 
sphere and its OH lines, and also by the optical faintness 
of the targets (/ > 23 mag). The high- redshift popu- 
lation of early-type galaxies is, thus, best investigated 
from space with low-resolution spectroscopy. This is the 
case of the HST Advanced Camera for Surveys (ACS), 
which is equipped with a slitless grism covering the spec- 
tral range 5500 - 10500 A with a dispersion of 40 A/pix 
in the first spectral order (Pasquali et al. 2006a). The 
HST/ ACS grism has been used for spectroscopic surveys 
like GRAPES (Pirzkal et al. 2004) and PEARS (Mal- 
hotra et al., in preparation), where early- type galaxies 
brighter than / ~ 26 mag have been identified up to 
z=2.5 and their individual spectra (typically with S/N > 
10 pixel -1 ) have been analyzed to constrain their stellar 
populations. The low resolution of these spectra clearly 
prevent us from measuring line indices, but allows us to 
detect the broad feature Mgw in the rest-frame 2640- 
2850 A (Daddi et al. 2005; Maraston et al. 2006) at 
z>1.3 and measure the 4000A break at 0.4<z<1.2 and 
use it along with the continuum to estimate stellar ages 
and mctallicities down to M* ~ 10 9 — 10 10 M Q (Pasquali 
et al. 2006b). This limiting stellar mass is a factor of 
about 10 lower than what can be achieved with ground- 
based spectroscopy, and makes it possible, for the first 
time, to study the stellar populations at the faint end 
of the red sequence at z~l. For example, Pasquali et 
al. (2006b) identified 18 early-types in the Hubble Ultra 
Deep Field (Beckwith et al. 2006) with 0.49 < z < 1.15. 
They fitted the individual HST/ ACS grism spectra with 
synthetic star formation histories, based on the models 
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Fig. 1. — Observed properties of the PEARS early- type galaxy 
sample, extracted from visual classification in GOODS N and S 
HST/ ACS images down to jab = 24 mag (full sample shown as 
dots, Ferreras et al. 2009a). We cross-correlate the total sample 
of 377 (CDFS) and 533 (HDFN) early-types with the PEARS ex- 
tractions. After visual rejection of bad spectra and rejection of 
redshfts z<0.4 we end up with the sample presented in this fig- 
ure (circles): 91+13 (CDFS) and 104+20 (HDFN) galaxies. The 
first number represents those with a "red" SED (grey filled cir- 
cles): prominent 4000-break and no emission lines. The second 
number corresponds to "blue" SEDs (open circles, see text for de- 
tails). We overlay an 'X' symbol on those galaxies with an X-ray 
detection from the Chandra deep fields (Alexander et al. 2003; 
Luo et al. 2008). Notice that all X-ray detections correspond to 
"red" galaxies. The dashed line in the bottom panel illustrates the 
resolution limit at 2R e ~ FWHMhst/ACS- The solid line is the 
surface brightness limit, considering the cosmological dimming, at 
«AB = 24 mag (AB), assuming a surface brightness at z=0.7 of 
MF775W = 22 mag/arcsec 2 . 

of Bruzual & Chariot (2003) and including chemical en- 
richment (Ferreras & Silk 2003). They found stellar ages 
varying from 3-4 Gyr at z>0.65 to 5-6 Gyr at lower red- 
shifts, stellar metallicities between 1 and 0.03 Z & and 
stellar masses between 3 x 10 9 and 3 x 10 n M Q . The 
estimated ages are consistent with a formation redshift 
z>2-5. 

In this paper, we extend the analysis of Pasquali et 
al. (2006b) to a sample of early-type galaxies that is 
a factor of 11 larger than in the HUDF. The objects - 
228 in total are extracted from the sample of mor- 
phologically selected early-type galaxies of Ferreras et 
al. (2009a) in the GOODS North and South fields and 
have been observed as part of PEARS 1 . This is the 
largest database of individual, high S/N grism spectra 
obtained so far for early- type galaxies at z~l, that al- 
lows a study of their stellar populations as a function of 

1 http : //archive . stsci . edu/prepds/pears/ 
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Fig. 2. — Rcst-framc properties of PEARS early-type galaxies. 
The sample is shown with respect to absolute V-band magnitude 
(top), stellar mass (middle) and projected half-light radius (bot- 
tom). The masses and magnitudes are computed using the best- fit 
SFH. 'Red' and 'blue' galaxies are represented by filled and open 
circles, respectively. The solid lines at the top panel track the 
«AB = 24 mag cut in apparent magnitude for two extreme cases: 
an old, passively evolving population (solid) or an extended star 
formation history (dashed). The small dots correspond to the orig- 
inal sample of visually-classified spheroidal galaxies from Ferreras 
et al. (2009a). We overlay an 'X' symbol on those galaxies with 
an X-ray detection from the Chandra deep fields North and South 
(Alexander et al. 2003; Luo et al. 2008). 

stellar mass and redshift, and can place more stringent 
constrains on the formation and evolution of (not only 
massive) early type galaxies since z~l. 

Throughout this paper a standard ACDM cosmology 
from the WMAP-year 5 data is used to determine ages 
and rest-frame properties (Komatsu et al. 2009). 

2. COLLECTING THE SAMPLE OF EARLY- TYPE 
GALAXIES 

The data consist of ACS images taken with the WFC 
G800L grism. The PEARS survey comprises eight point- 
ings towards the GOODS North and South fields - 4 
pointings each - with three HST roll angles for each 
pointing (20 HST orbits per field), plus the HUDF with 
four roll angles (40 HST orbits). A forthcoming paper 
(Malhotra et al. 2009, in preparation) will describe the 
PEARS project and data in detail. The HUDF data 
comes from an earlier similar project, GRAPES (Pirzkal 
et al. 2004; Pasquali et al. 2006b). We use the PEARS 
IDs from the master catalogue of source extractions. The 
G800L grism produces slitless, low-resolution spectra in 
the wavelength range 6000-9500A. The optimal spectral 
resolution (R ~ 100) is achieved for unresolved sources, 
whereas for our sample of early-type galaxies, the effec- 
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TABLE 1 

Catalog of PEARS early- type galaxies: Observed data 
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RA 


Dec 
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0.49 


R 


52017 
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+62.208328 


23.35 


3.78 


0.97 


R 


61235 


53.1881218 


-27.827766 


22.87 


6.35 


1.06 


R 


69419 


53.1323776 


-27.814238 


22.91 


5.72 


0.75 


R 


40498 


189.1528931 


+62.186619 


23.72 


7.77 


0.85 


B 


47252 


189.2056274 


+62.198719 


23.41 


3.60 


0.96 


B 


42729 


53.2094307 


-27.861725 


23.79 


5.44 


0.56 


B 


17362 


53.1624603 


-27.914028 


22.49 


6.72 


0.51 


B 


57554 


189.2057190 


+62.219887 


22.21 


5.83 


0.48 


B 


119723 


189.3504181 


+62.327610 


22.89 


4.12 


0.80 


B 


53462 


53.1896324 


-27.841940 


22.79 


5.84 


0.56 


B 


107477 
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-27.739756 


22.29 


13.32 


0.42 


B 


104645 
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+62.301193 


23.78 
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B 



Note. — A full version of tabic 1 is published in its entirety in the 
electronic edition of the Astrophysical Journal. A portion is shown here 
for guidance regarding its form and content. 

a PEARS ID number. 

k'Red' vs 'Blue' early-type (see text for details). 
c l pixel — 0.03 arcscc 



tive resolution is lower, around R ~ 50. 

The master catalogue of PEARS detections was cross- 
correlated with the sample of 910 visually-classified 
early-type galaxies in the GOODS North and South fields 
of Ferreras et al. (2009a). This sample completes the 
one presented previously in the CDFS field (Ferreras et 
al. 2005). The selection consists of visual inspection 
(in all four ACS bands) of all source detections brighter 
than iF775W < 24 mag (AB) by four people, with two 
rounds of visual inspection for those galaxies that were 
harder to classify. No colour cuts or CAS-based selec- 
tions were made, although the sample is consistent with 
the high concentration and low asymmetry expected in 
spheroidal galaxies with respect to later morphological 
types (Lisker et al., in preparation). Furthermore, the 
visually-classified sample comprises 80% of red, quies- 
cent stellar populations. We believe visual inspection 
alone is much more efficient at selecting spheroidal galax- 
ies, eliminating contamination from late-type galaxies or 
unresolved sources, and removing the bias inherent in 
color-based selection methods. From this catalogue we 
extract the available PEARS spectra of each identified 
object from our database. The spectra are background 
corrected and scaled to match the direct image magni- 
tudes of each object. Each file contains the spectra ex- 
tracted for the different PAs (typically around three PAs 
per galaxy). We compute a simple average of the spec- 
tra available for each object, possibly excluding those 
PAs whose spectra are clearly deviant. The errors in the 
final combined spectra are computed by propagating the 
flux errors in each individual SED. The spectra were vi- 
sually inspected, and some of them had to be discarded 
because of a very low signal-to-noise ratio, contamina- 
tion from nearby sources, or because they were located 
close to the edge of the field. Out of the original sample 
of 910 galaxies (533 in HDFN and 377 in CDFS) we ex- 



tracted 136 spectra in the North and 131 in the South. 
The number of rejected galaxies is 19 and 18 in the North 
and South fields, respectively. 

Regarding rcdshift, we start with estimates from the 
PEARS-based general catalog of photometric redshifts 
(Cohen et al. in preparation), and the photometric red- 
shift catalogue of Mobasher et al. (2004) and Ferreras 
et al. (2009a). Spectroscopic redshifts are available 
from FORS2 (Vanzella et al. 2008) and the VIRMOS- 
VLT Deep Survey (VVDS, Popesso et al. 2008) in the 
South, and from the Team Keck Redshift Survey (TKRS, 
Wirth et al. 2004) in the North. However, our slit- 
less grism data has enough S/N to allow for an estimate 
of the redshifts directly from the low-resolution spectra. 
Two of us (I.F. and A. P.) inspected visually each SED 
and compared it with a number of template spectra to 
determine a rough estimate of the rcdshift. Then we 
ran a simple code that uses the guess and templates to 
accurately determine the redshift via a standard maxi- 
mum likelihood method. The accuracy of our redshifts 
is as good as spectroscopic estimates for most of the 
galaxies, especially those with a prominent 4000 A break 
(i.e. the majority of the sample) or those with emis- 
sion lines. Incidentally, some of our (grism-based) red- 
shift estimates improved on the spectroscopic values from 
pipeline-processed databases such as GOODS/FORS2. 
This improvement is to be expected for the faintest galax- 
ies with a quiescent spectra. In those cases the contin- 
uum from the ground-based spectroscopic data is very 
noisy and, simply put, our low-resolution spectra yield 
more photons per pixel over the CCD, significantly im- 
proving the signal. Hence, we consider our quoted red- 
shifts as spectroscopically accurate. 

Finally, we apply a cut in redshift, removing from the 
analysis all galaxies with rcdshift z<0.4. At those red- 
shifts the 4000A break - on which the modelling strongly 
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Fig. 4. — Rcst-framc U — V color-magnitude relation (right) and 
color-stellar mass relation (left). Filled (open) circles correspond 
to photometric 'red' ('blue') galaxies. The dashed line gives the 
z~ colour-magnitude relation of Coma early-type galaxies from 
Bower et al. (1992). 



Fig. 3. — Comparison of half-light radius (R e ; top) and V- 
band M/L (Ty with respect to the solar value; bottom) with a 
z~ sample of galaxies with Sersic index ng > 2.5, extracted 
from the SDSS NYU Value Added Galaxy Catalogue (grey dots, 
Blanton et al. 2005). The half-light radius (top) and M/L in the 
V band is shown as a function of stellar mass (left) or absolute 
magnitude. 'Red' ('blue') PEARS galaxies are represented as black 
solid (open) circles. 

relies - falls outside of the sensitivity region of the G800L 
grism. This cut reduces the sample to 124 and 104 galax- 
ies in the North and South fields, bringing the total to 
228 early-type galaxies with z>0.4 and ^F775W < 24 mag 
(AB). The median redshift of the sample is zm=0.75. If 
we measure the signal-to-noise ratio as the average S/N 
per pixel within the spectral window 8000-8500A , our 
sample has a median value of S/N~ 15 per pixel. 

2.1. Sample Properties 

The main observational data are presented in table 1. 
Figure 1 shows the distribution of total apparent magni- 
tude (top) and size (half-light radius; bottom). The orig- 
inal sample of visually-classified galaxies from Fcrreras et 
al. (2009a) is shown as small dots. The PEARS sample 
of early-type galaxies is shown as circles. We subdivide 
the sample into 'red' (grey dots) and 'blue' types (open 
circles). This classification is taken from the original data 
and is determined by the spectral template that gives the 
best fit in the analysis of photometric redshifts. For the 
PEARS sample presented here, the division is much more 
accurate because this information is complemented by 
the visual appearance of the low-resolution SED. 'Red' 
galaxies feature prominent 4000A breaks, whereas 'blue' 
galaxies have a clear blue continuum with the occassional 
presence of emission lines. The sample comprises 195 
'red' and 33 'blue' spheroidal galaxies. The accuracy 
of the total magnitudes and half-light radii is estimated 
from simulations of synthetic galaxies with Sersic profiles 
with similar values of size and magnitude as the original 
sample. We recover the original values with an uncer- 
tainty Ai F 775w = 0.05 mag and AR e /R e = 0.09 (Fer- 
reras et al. 2009a). We overlay 'X' symbols over those 
galaxies with an X-ray detection from the Chandra Deep 
fields North and South (Alexander et al. 2003; Luo et 
al. 2008). 22 sources have an X-ray detection which ex- 
tends over the full range in apparent magnitude. It is 
worth emphasizing that all X-ray detected sources in our 
sample are photometrically classified as 'red' galaxies. 

The distribution of rest-frame properties of our sam- 



ple is shown in figure 2, where the absolute magnitude 
in the V band (top); stellar mass (middle) and physical 
half-light radius (bottom) is shown with respect to red- 
shift. The K-correction and the M/L estimate needed 
for the top two panels is taken from the best fit models 
of the runs explained in §3. The 'red' ('blue') galaxies 
are shown as filled grey (open black) circles. 'Blue' early- 
types are mostly faint, low-mass galaxies. We should em- 
phasize that our selection process (visual classification) 
does not introduce any bias regarding the color distri- 
bution. The lines in the top panel correspond to the 
^F775W < 24 mag limit imposed on the original sample. 
Since the translation from apparent to absolute magni- 
tude depend on the template SED for the K-correction, 
we show the apparent magnitude limit corresponding 
to two types of populations: quiescent (solid line) and 
young (dashed line) . These two cases are the extrema of 
the photometric types used in Ferreras et al. (2009a), 
namely t=0 and t=7. These types correspond to an ex- 
ponential star formation history at fixed solar metallicity, 
beginning at a formation redshift zp = 3 with timescale 
r = 0.05 Gyr (t=0) or r = 8 Gyr (t=7). Similarly to fig- 
ure 1, the 'X' symbols are overlaid on those galaxies with 
an X-ray detection from the Chandra Deep Fields North 
and South (Alexander et al. 2003; Luo et al. 2008). Al- 
most all X-ray detections correspond to galaxies brighter 
than My < —20, with stellar masses above 10 10 M Q . 
However, many of them are concentrated, with half-light 
radii R e ~ 1 kpc. 

Figure 3 shows the physical half-light radius (R e ; top) 
and M/L in the V band (Ty with respect to the so- 
lar value; bottom) as a function of stellar mass (left) or 
absolute magnitude (right). Black filled (open) circles 
correspond to 'red' ('blue') galaxies. In order to com- 
pare with z=0 galaxies, we include as grey dots the val- 
ues from the DR4 NYU Value Added galaxy catalogue 
(Blanton et al. 2005) with a Sersic index ns > 2.5. Local 
galaxies appear significanty larger and with higher M/L. 
The latter can be explained by pure passive evolution of 
the stellar populations, but the change in size requires 
a mechanism beyond stellar evolution. This result is 
shown for the full sample of GOODS spheroidal galaxies 
in Ferreras et al. (2009b) and confirms previous results 
at similar or higher redshifts (Trujillo et al. 2007). Size 
evolution is even stronger when galaxies at z~2 are con- 
sidered (van Dokkum et al. 2008; Buitrago et al. 2008). 
Recent attempts at explaining the size evolution invoke 
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Fig. 5. — Quantitative morphology of the PEARS early-type 
galaxy sample. The plot shows the concentration, M20 (a normal- 
ized second order moment of the 20% brightest pixels) and asym- 
metry. Filled grey (black open) circles correspond to 'red' ('blue') 
spheroidal galaxies. See text for details. 

a significant ejection of gas mass caused by quasar feed- 
back (Fan et al. 2008), or via dry mergers in the standard 
paradigm of hierarchical evolution (Hopkins et al. 2009). 

The rest- frame U — V color- magnitude and color-stellar 
mass relations are presented in figure 4, with the same 
coding for photometric type as in the previous figures. 
Notice color is correlated more strongly with stellar mass 
than with absolute magnitude. The dashed line gives the 
Coma color-magnitude relation from Bower et al. (1992) 
for E+SOs. The extent of the line covers the magnitude 
range probed in Coma. The relation at high redshift is in 
good agreement with Coma if we consider that passive 
evolution of a stellar population with solar metallicity, 
formed at zf = 3 changes its U— V colour by ~ 0.3 mag 
between the median redshift of our sample (z) <~ 0.7 and 
z= 0. Instead of dealing with simple, K-corrected colors, 
we use the spectroscopic data to determine the best fit 
ages from models of stellar populations. In the next sec- 
tion we compare the low-resolution SEDs from PEARS 
with a large grid of star formation histories to transform 
this information into a more meaningful diagram involv- 
ing stellar ages. 

A rough morphological quantification of the sample 
is presented in figure 5, where the concentration, M 2 o 
and asymmetry values are shown for the zf775W images 
(Lisker et al. in preparation). Concentration (see e.g. 
Bershady et al. 2000) is defined as C= 5 logRgo/R^o , 
where Kx is defined as the radius within which the en- 
closed flux corresponds to X % of the total flux (as mea- 
sured within the Petrosian radius). M 2 o is defined as the 
normalized second order moment of the brightest 20% of 
the pixels no matter where they are (Lotz et al. 2004). 
The asymmetry parameter (Aigo) is obtained by sub- 
tracting a 180° rotated image from the original. Ai 80 
is given by this difference - summed and expressed as 
a fraction of the original flux (Conselice 2003). The fig- 
ure shows that most galaxies (whose selection was purely 
based on a visual classification, i.e. NOT extracted from 
their CAS values) are highly concentrated both in C and 
M 2 o and have low asymmetry. Notice the blue photomet- 
ric types (open circles) coincide with the lowest values 
of concentration, although those values are still high for 
a general population involving all galaxy morphologies, 
which can have concentrations as low as C~1.5. The 



TABLE 2 

Range of parameters explored for the three models 
considered in this paper (see text for details). 



Model/Param 


MIN 


MAX 


Comments 


SSP 

t 

log(Z/Z ) 


0.1 Gyr 
-1.5 


2 params 
z F = 10" 
+0.3 


Age 
Metallicity 


EXP 

logr (Gyr) 
Z F 

log(Z/Z ) 


-1 

z-0.1 Gyr 6 
-1.5 


+0.6 

10 
+0.3 


3 params 
Exp. Timescale 
Formation epoch 
Metallicity 


2BST 

to 
ty 
fy 

log(Z/Z ) 


2 

0.1 
0.0 
-1.5 


z F = 10" 
2 

1.0 
+0.3 


4 params 
Old (Gyr) 
Young (Gyr) 
Mass fraction 
Metallicity 


CSP 

v (Gyr" 1 ) 
Tf (Gyr) 
$ (Gyr) 

ZF 


20 

-1 


z-0.1 Gyr a 


+0.3 
1 

10 


3 (free) params 
SF Efficiency (fixed) 
Gas Infall timescale 
Outflow fraction 
Formation epoch 



a SSP,2BST: The oldest age possible for the old component corresponds 
to a formation redshift zp — 10 at the redshift of the galaxy, e.g. 6.7 Gyr 
for a z— 0.7 galaxy. 

kEXP,CSP: The latest formation redshift available corresponds to the 
observed redshift of the galaxy minus 0.1 Gyr, e.g. zf > 0.72 for a z— 0.7 
galaxy. 



asymmetry of blue and red early- types is similar. 

3. MODELLING THE STAR FORMATION HISTORIES OF 
GALAXIES 

Our slitless spectra have good enough S/N down to the 
*F775W = 24 mag (AB) limit of the original sample. The 
S/N ranges from 30 down to 4-5 per pixel at the limit- 
ing magnitude. One of the main advantages of G800L 
grism data for the analysis presented here is the superb 
flux calibration of the spectra enabled by the optimal 
flat fielding of the images. We estimate flux calibration 
systematic errors below 5% within the spectral range of 
interest (Pirzkal et al. 2004). 

The wavelength coverage of the ACS/G800L grism 
(6000-9500A) is ideally suited for the analysis of old 
stellar populations at z=0.4-1.5, because the prominent 
4000A break falls within its sensitivity range. Further- 
more, the compact nature of early-type galaxies and their 
homogeneous stellar populations minimises the contami- 
nation and loss of spectral resolution compared to slitless 
grism data of late-type or irregular galaxies. 

3.1. Methodology 

We follow a standard likelihood method fitting the low- 
resolution spectra to grids of models. A proper extrac- 
tion of the age and metallicity distribution requires the 
use of composite models of stellar populations. Simple 
stellar populations (i.e. models with a single age and 
metallicity) will give luminosity-weighted age estimates 
that can differ significantly from a more physical mass- 
weighted age according to a composite stellar population 
(Ferreras & Yi 2004; Pasquali et al. 2005; Serra & Tragcr 
2007; Rogers et al. 2008). In order to assess systematic 
differences related to the parameterisation of the star for- 
mation histories, we consider the following four models 
(whose parameters are summarised in table 2). 
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Fig. 6. — A sample of best fits to the grism data of 'red' early- 
type galaxies using models with chemical enrichment (CSP). Fits 
to the other two models considered in this paper (EXP and 2BST) 
do not differ much from these. The observed SEDs are shown in 
grey, including a la error envelopes (dashed lines). The model fits 
are shown as black solid lines. Each galaxy is labelled (from top 
to bottom) by their PEARS ID, iF775w magnitude and redshift. 

I. Simple Stellar Populations (SSP): SSPs are the 
building blocks of any population synthesis model. 
They correspond to a stellar population with a sin- 
gle age and metallicity. Although SSPs are often a 
fair approximation of a globular cluster, the longer 
timescales and more complex chemical enrichment 
history of a galaxy makes an SSP a very rough 
approximation in this case. Furthermore, the fact 
that a single age and metallicity are used to fit 
the data can introduce a significant difference be- 
tween the SSP-determine age and a more physical 
mass-weighted age. Nevertheless, we include these 
models in the analysis for comparison, given that 
they are the simplest models from which the com- 
posite models are generated. This case explores a 
grid of 128 x 128 SSPs over a wide range of ages 
and metallicities, as given in table 2. 

II. r-model (EXP): The SFH is described by an ex- 
ponentially decaying function, with a timescale (r) 
and a formation redshift (zp), corresponding to the 
epoch at which star formation starts. The metal- 
licity is the third free parameter and it is kept con- 
stant throughout each SFH. 

III. 2-Burst (2BST): Superposition of two simple stellar 
populations, described by four parameters: the age 
of the old (to) and the young components (ty), the 
mass fraction in young stars (fy) and the metallic- 
ity of the system (assumed to be the same for both 
populations). 

IV. Chemical Enrichment Model (CSP): We follow a 
consistent chemical enrichment code that incorpo- 
rates the evolution of metallicity with the star for- 



FlG. 7. — Same as figure 6 for a sample of 'blue' early-type 
galaxies. The fits also correspond to CSP models. In this case the 
fitting method excludes spectral windows of size ±100A around 
[O III], [O II], H/3, and Ho lines whenever the redshift moves 
them into the spectral range of the grism data. (See Straughn ct 
al. (2008) for details of emission line galaxies in the PEARS HUDF 
spectra). 

mation rate. The model is defined in Ferreras & 
Silk (2000) and was previously applied to similar 
data from the GRAPES sample of early-types and 
bulges in the HUDF (Pasquali et al. 2006b; Hathi 
et al. 2009). The formation history is described 
by four parameters: a star formation efficiency (v ) 
controlling the transformation of gas into stars via 
a Schmidt law; an outflow parameter ((3) which 
is the mass fraction of gas ejected i\ir> <X result of 
the ongoing star formation rate; the exponential 
timescale of gas infall (r/); and the formation epoch 
given by the redshift (zp) at which infall starts. 

For a given choice of parameters, the models determine 
the distribution of stellar ages and metallicities. These 
distributions are used to combine simple stellar popula- 
tions from the latest 2007 models of Chariot & Bruzual 
(hereafter CB07, Bruzual & Chariot 2003; Bruzual 2007, 
latest models in preparation) , assuming a Chabrier initial 
mass function (Chabrier 2003). Notice the CB07 mod- 
els have been recently updated with respect to the high 
resolution spectra. However, at the low resolution grism 
spectra considered in this paper, the difference between 
CB07 and the expected CB09 models is negligible (Char- 
lot, priv. comm.). The resulting spectra are degraded to 
the slitless PEARS SED, and compared via a standard 
maximum likelihood method. Each model is explored on 
a grid comprising between 32,000 and 65,000 realizations 
for a range of parameters as described in table 2. The 
best fit model is used to determine the age and metallic- 
ity of the underlying populations, and the M/L required 
to translate flux into stellar mass. 

3.2. Extracting star formation histories 
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Figures 6 and 7 show some of the best fits of the CSP 
models for the 'red' and 'blue' galaxies, respectively. The 
model SED is shown as a solid black line, and the ob- 
served SED is shown in grey (including the la uncer- 
tainty as dashed lines). Notice 'blue' galaxies have emis- 
sion lines, most notably [O II], [O III] and H/3. These 
lines are masked out in the likelihood analysis. The 
blue spectra show a significant range of ages, from large 
4000A breaks (as in PID 119723) to very young stars (e.g. 
47252). In this paper the analysis only deals with stellar 
populations, hence the presence of an AGN could sig- 
nificantly alter the estimated ages of the 'blue' galaxies. 
Nevertheless, figure 2 shows that the subsample of 'blue' 
galaxies is dominated by low-mass systems, for which 
one would expect a significant contribution to blue light 
from young stellar populations. One could argue that 
these low stellar mass estimates are biased because of the 
contribution from the AGN. However, the apparent mag- 
nitude of all 'blue' galaxies in our sample is consistently 
fainter than 'red' galaxies (see figure 1). The 'red' sub- 
sample (figure 6) looks much more homogeneous, with a 
prominent 4000Abreak and no emission lines. The SEDs 
also display the characteristic dip around the G-band 
(rest-frame 4300A) and for the lower redshift galaxies one 
can also discern the Mg feature at rest-frame 5170A(see 
e.g. PID 122743 or 127697 in figure 6). Needless to 
say, the low resolution of the G800L grism prevents us 
from doing any analysis of the absorption lines. How- 
ever, the "smoothed-out" features of these spectra can 
contribute to constraining the age and metallicity dis- 
tribution. Indeed, higher resolution spectra (R/^2000) 
of Virgo elliptical galaxies at very high signal-to-noise 
ratio (S/N> lOOA -1 ) has shown that independent fits 
to the equivalent widths or a direct fit to the SED give 
age distributions that are in good agreement (Rogers et 
al. 2008). 

4. RESULTS AND CONCLUSIONS 

Figure 8 shows the (mass- weighted) average age (top), 
age scatter (defined as an RMS, middle) and metallic- 
ity (bottom) as a function of stellar mass for three bins 
in redshift (top panel). These estimates correspond to 
the chemical enrichment models (CSP). The best fit star 
formation history from the models is used to compute 
the average and RMS of the age distribution. 'Red' 
and 'blue' early-types are shown as filled grey or open 
black circles, respectively. Characteristic error bars are 
included. The figure includes simple fits to the age dis- 
tribution for the lowest redshift sample (black line in the 
top-left panel) and the highest redshift bin (black line in 
the top-right panel). The lines in grey are just copies of 
these two lines in the other panels, to guide the eye. The 
arrow in the top-right panel is the expected age differ- 
ence between the high and the low redshift bin (i.e. pure 
passive evolution). 

The figure allows us to understand the blue cloud- 
red sequence diagram of the sample. Our blue-cloud 
spheroidal galaxies are younger systems with roughly 
similar metallicities to the red sequence distribution. It 
is worth mentioning that the CSP models suggest all our 
galaxies have a similar spread of stellar ages, and it is 
only the average age that correlates with stellar mass. 
This result is consistent with the more detailed analysis 
of Rogers et al. (2008) on high-S/N spectra of cllip- 
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Fig. 8. — Best fit average age (top), age spread (middle, defined 
as the RMS of the distribution) and metallicity (bottom) for CSP 
models. The lines in the top panel are simple least square fits to the 
red galaxies for the lowest and highest redshift bins. These lines 
are shown in black for the subsamples where the fit was done. The 
grey lines are copies of these two fits, to guide the eye. The arrow 
in the top right panel corresponds to the lookback time between 
the median redshift of the high- and low redshift bins. On the left 
panels, characteristic error bars are shown. 'Red' ('blue') early- 
type galaxies are shown as grey filled (black open) circles. 



tical galaxies in the Virgo cluster. The lines shown in 
the top panel illustrate that passive evolution over this 
redshift range can explain the 'red' spheroidals - which 
represent the large majority of the sample. The average 
ages of the subsample at lower redshift (0.35 <z< 0.67; 
top-left) suggest that the correlation between stellar age 
and mass is such that instead of a monotonic trend, the 
scatter in age increases in galaxies with stellar masses 
below 10 10 M (see e.g. Caldwell et al. 2003; Gallazzi 
et al. 2006; Yamada et al. 2006; Rogers et al. 2008). 

In order to check systematic effects arising from the 
specific way to parameterise the models, we compare in 
figure 9 the average age and metallicity distribution of 
the three models considered. There is overall good agree- 
ment between the CSP and the EXP models (top pan- 
els), whereas the age and metallicities of SSP and 2BST 
models differ significantly with respect to CSP models. 
Notice that best fit SSP metallicities are often very low, 
as the best fit wanders along the age-metallicity degen- 
eracy, affecting the estimated ages. Composite models 
have been found to give more realistic metallicity con- 
straints than SSPs (Ferreras & Yi 2004). However, we 
cannot look into the reduced x 2 values to rule out one 
model against another: all models give equally accept- 
able values of Xr ~ 1- Observations at higher spectral 
resolution dot not fare better at disentangling this de- 
generacy (Rogers et al. 2008). 

It is also worth mentioning that for the composite 
models (2BST, EXP and CSP) there is good agreement 
between the best-fit metallicities for the red galaxies - 
within a realistic ~ 0.3 dex uncertainty. Blue galaxies, as 
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TABLE 3 

Catalog of PEARS early-type galaxies: Physical properties (CSP models) 



PID 


M s 


U-V 


M v 


R e 


(t) 


of 


(log ZlZn,) 


8 




log T 




x 1 10 M^ 






kpc 


Ovr 


Ovr 








Ovr 


65620 


11.51 


1.06 


-23.20 


4.13 


3.82 


0.81 


-0.18 


0.28 


7.5 


-0.48 


56798 


11.02 


1.28 


-21.17 


1.98 


6.21 


1.12 


-0.10 


0.21 


9.8 


-0.31 


83499 


11.24 


1.05 


-22.34 


3.06 


4.40 


0.81 


-0.28 


0.37 


7.5 


-0.56 


75762 


11.26 


1.24 


-22.43 


3.56 


4.73 


0.92 


-0.15 


0.24 


9.9 


-0.35 


61447 


10.76 


1.14 


-20.83 


0.95 


4.39 


0.86 


-0.20 


0.32 


7.6 


-0.47 


127697 


10.29 


1.17 


-19.39 


1.31 


5.81 


1.06 


-0.11 


0.17 


5.5 


-0.53 


52017 


10.19 


0.80 


-20.95 


0.90 


2.71 


0.71 


-0.30 


0.41 


3.2 


-0.50 


61235 


10.86 


1.06 


-22.16 


1.55 


3.64 


0.77 


-0.22 


0.33 


9.8 


-0.48 


69419 


10.35 


1.08 


-20.42 


1.26 


4.51 


0.89 


-0.18 


0.26 


6.0 


-0.48 


40498 


9.33 


0.02 


-19.47 


1.79 


1.05 


0.53 


-0.36 


0.45 


0.8 


-0.07 


47252 


8.69 


-0.49 


-19.87 


0.86 


0.10 


0.04 


-0.94 


0.41 


0.9 


-0.26 


42729 


8.83 


0.37 


-18.32 


1.06 


0.66 


0.38 


-0.38 


0.45 


0.7 


-0.12 


17362 


9.72 


0.47 


-19.45 


1.24 


2.35 


0.74 


-0.47 


0.58 


1.2 


-0.22 


57554 


9.85 


0.48 


-19.58 


1.04 


1.98 


0.71 


-0.54 


0.61 


1.0 


-0.19 


119723 


9.65 


0.09 


-20.15 


0.93 


0.65 


0.45 


-0.13 


0.20 


1.1 


0.10 


53462 


9.45 


0.26 


-19.36 


1.13 


0.92 


0.50 


-0.40 


0.48 


1.0 


-0.04 


107477 


10.31 


1.41 


-19.16 


2.21 


2.27 


0.75 


-0.36 


0.46 


1.0 


-0.22 


104645 


8.79 


0.38 


-17.93 


1.00 


0.95 


0.51 


-0.35 


0.44 


0.6 


-0.07 



Note. — A full version of tabic 3 is published in its entirety in the electronic edition of the Astrophysical 
Journal. A portion is shown here for guidance regarding its form and content. 

a 0-t is the root mean square of the stellar age distribution. 




Fig. 9. — Comparison between the predicted average ages and 
mctallicities for the models explored in this paper (see text for 
details). The solid line gives the 1:1 mapping. Solid (open) circles 
correspond to 'red' ('blue') early-type galaxies. 

expected, give more discrepant results, mainly caused by 
the lack of spectral features in the low-resolution PEARS 
SEDs which makes metallicity estimates more uncertain. 

Another check of our results was done with respect to 
the older version of the population synthesis models of 
Bruzual & Chariot (2003). In figure 10 we compare the 
average age and metallicity when using either the 2003 
models (horizontal axes) or the improved CB07 models 
(vertical axes). The discrepancy between both versions 
is most prominent for the 'blue' subsample (open cir- 
cles), which reflects the different treatment of some of 



Fig. 10. — We illustrate here the effect of a different choice of 
stellar population models. The average age (top) and metallic- 
ity (bottom) of our models, based on the latest Chariot & Bruzual 
(2007) models, is compared with the predictions when the previous 
version of population synthesis models arc used (Bruzual & Char- 
lot 2003). Notice 'red' galaxies (filled circles) give similar results, 
whereas the younger populations found in the 'blue' subsample 
(open circles) depend strongly on the choice of population models. 

the young-intermediate age phases of stellar evolution. 
The difference of age or metallicity between BC03 and 
CB07 models for older populations ('red' galaxies, grey 
circles) is negligible within the error bars. 

The average and spread of the age distribution shown 
in figure 8 can be shown with respect to the phenomeno- 
logical parameters used to describe the star formation 
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Fig. 11. — Best-fit parameters for a continuous model including 
chemical enrichment (CSP, top) or a two-burst scenario (2BST, 
bottom). The CSP parameters are (from left to right), the outflow 
fraction (f3), formation epoch (zp; defined as the time when gas 
infall starts) and the exponential timescale of gas infall (tj ). The 
2BST parameters are (also from left to right) the mass fraction in 
young stars (/y); young stellar age (*y) ; and old stellar age (to)- 



history. Figure 11 shows the best-fit parameters for the 
CSP (top) and the 2BST (bottom) models as a function 
of stellar mass. A large spread is found for the fraction of 
gas and metals ejected in outflows (/?), whereas the gas 
infall timescale is overall rather short (tj < 1 Gyr). Over 
the 0.4<z<1.2 range probed by our sample - correspond- 
ing to a cosmic age between 5 and 10 Gyr (with 13.5 Gyr 
being our present age) - it is the formation epoch (top 
middle panel) what mainly correlates with stellar mass 
to give the age- mass trend shown in figure 8. Our 'blue' 
early- types have formation epochs zp ~ 1, which - given 
the redshift of the sample - imply a large amount of 
young stars. 

If we consider instead a two-burst model (2BST; bot- 
tom panels), the main driver of the age-mass relation is 
given by the young stellar mass fraction (iy), although 
the mass fraction in young stars (/y) also shows a corre- 
lation with mass. 'Blue' early-types have high values of 
/y. In general, we find that significant masses in young 
stars (open circles) only appear in early-type galaxies 
with masses < 10 1 M©. However, notice there is a frac- 
tion of 'red' early-type galaxies - which live on the red 
sequence (figure 4) - with a significant amount of young 
stars (/y > 0.4). This is consistent with the observed in- 
crease of the recent star formation in spheroidal galaxies 
with redshift (Kaviraj et al. 2008). 

Both a smooth description of the star formation history 
(CSP) and a more abrupt model (2BST) confirm that 
red-sequence spheroidal galaxies with masses above ~ 
10 n M© have the bulk of their stellar populations formed 



at zp > 3. This result confirms previous findings (see 
e.g. Tantalo et al. 1996; Vazdekis et al. 1997; Perez- 
Gonzalez et al. 2008; Wiklind et al. 2008; Mancini et 



al. 2009) although our work improves on the analysis, 
since we use SEDs instead of photometry to constrain 
the age distribution, and we consider a large volume of 
parameter space to describe the star formation history. 

The old populations found in massive early-type galax- 
ies in combination with the lack of number density evo- 
lution (see e.g. Ferreras et al. 2009b) confirms that new 
channels must be involved to explain the presence of 
galaxies on the red sequence at the massive end. In anal- 
ogy to the cartoon version of the galaxy evolution tracks 
A-C in figure 10 of Faber et al. (2007), we propose a 
new one (track D) in the right panel of our figure 12. 
Tracks A-C (left panel) follow the standard growth via 
wet mergers on the blue cloud, followed by quenching of 
star formation (the near vertical black arrow) after which 
growth to the top of the red sequence consists of dry 
mergers (white arrow). This option is still controversial, 
since different semi-analytic models predict a very dif- 
ferent redshift evolution of the number density (Ferreras 
et al. 2009b, Hopkins et al. , in preparation). On the 
other hand, the recent work of Dekel et al. (2009) based 
on the Mare Nostrum cosmological simulations suggest 
a significant amount of massive galaxy growth via cold 
accretion at the intersection of cosmic filaments. This 
process would bypass the evolution in the blue cloud, as 
visualised by our track D on the right panel of figure 12. 
In this case, most of the stellar component in a massive 
galaxy would grow in situ, over a short and early period, 
explaining the old populations found in massive galax- 
ies (i.e. the "explosion" symbol in our track D). Fur- 
thermore, this process would be very efficient, turning 
gas into stars over short timescales, thereby explaining 
the small infall timescales found (figure 11). The fact 
that formation epoch correlates with stellar mass would 
indicate the connection between the stellar mass of the 
galaxy and the significance of the overdensity from which 
these galaxies grow. Lower masses correspond to 'lower- 
a' fluctuations thereby collapsing later. Nevertheless, 
it is worth emphasizing that the channel proposed by 
Dekel et al. (2009) would dominate at the most massive 
end. Intermediate mass galaxies will have a more com- 
plex history, with both tracks A-C and D contributing to 
their formation. The picture is unfortunately not fully 
explained, as massive galaxies also experience a very sig- 
nificant size evolution between z~2 and z=0 (Daddi et 
al. 2005; van Dokkum et al. 2008; Buitrago et al. 2008), 
an issue not explained yet in the context of Dekel et 
al. (2009). 

As a final note, we may speculate on the importance of 
the formation redshifts shown in the top-middle panel of 
figure 11. The best-fit values range between zp ~0.5 to 
zp ~10, as expected in hierarchical formation scenarios 
that produce a broad peak in the cosmic SFH at z~l-2. 
However, we also note that the highest formation red- 
shifts pile up either between zp ~7-8 or around zp ~10. 
While the estimates of the current spectrophotometry 
and SED modelling is undoubtedly not accurate enough 
to state the difference between the epochs zp ~7-8 and 
zp ~10 with any great confidence, it is tempting to no- 
tice that these redshifts coincide remarkably well with 
the epochs in WMAP-year5 cosmology where the peak 
of the Pop III star-formation should have occurred at 
z~10.8±1.8 in order to produce the WMAP-year5 polar- 
ization optical depth (Dunkley et al. 2009), and with the 
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Fig. 12. — Our proposed additional channel for the formation of galaxies on the red sequence is shown as 'D track' (right), in comparison 
with the standard picture from Fabcr et al. (2007) (A-C track, left). The age distribution of our sample of moderate redshift early-type 
galaxies along with recent findings of a lack of evolution of the comoving number density or the intrinsic colour distribution of massive 
galaxies (see e.g. Fontana et al. 2006; Conselice et al. 2007; Ferreras ct al. 2009b) suggest an in situ build-up of the stellar populations 
bypassing any trajectory in the blue cloud (except for the very short-lived formation stage corresponding to an intense, highly efficient 
process of star formation). This proposed channel can be justified by the cold accretion streams found in recent numerical simulations 
(Dekel et al. 2009). 



epoch of zp ~7-8 in which the global component of Pop II 
stars should have started shining in dwarf galaxies in or- 
der to finish reionization by z~6. Future data and mod- 
eling should concentrate on whether these epochs can be 
defined any more accurately than possible here. This 
analysis may be feasible with the upcoming HST/WFC3 
grism spectra of the earliest type galaxies at l<z<3 in 
the near future. 
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